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We use simulations to examine hysteresis and noise in a model system that produces heterogeneous 
orderings including stripe and clump phases. In the presence of a disordered substrate, these hetero- 
geneous phases exhibit 1/ f noise and hysteresis in transport. The noise fluctuations are maximal 
in the heterogeneous phases, while in the uniform phases the hysteresis vanishes and both a and the 
noise power decrease. We compare our results to recent experiments exhibiting noise and hysteresis 
in high-temperature superconductors where charge heterogeneities may occur. 

PACS numbers: 74.25.-q, 72.70.+m, 71.45.Lr 
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There is growing evidence that a wide variety of con- 
densed matter systems intrinsically exhibit heteroge- 
neous charge ordering in the form of mesoscopic clump, 
labyrinth, or stripe phases. Examples of systems where 
such phenomena may occur include cuprate supercon- 
ductors [1-5], antiferromagnetic insulators [6], and two- 
dimensional electron gas (2DEG) systems [7,8]. Hetero- 
geneous states are observed in manganites and diluted 
magnetic semiconductors [9] as well as superconductors 
[11,12]. In most of these systems, some form of quenched 
disorder from intrinsic defects in the sample is present, 
and it can destroy any long range ordering in the pat- 
terns [10]. Charged stripe forming systems have also 
been shown to exhibit self-generating disordered glassy 
properties [13-15]. The presence of ordered or disordered 
heterogeneities should affect the bulk transport, fluctu- 
ations, and transient responses of the systems; however, 
little is known about how the effect of heterogeneity on 
transport would differ from that of homogeneous but dis- 
ordered systems which form uniform crystalline or par- 
tially crystalline phases. 

Recently, transport experiments in underdoped sam- 
ples of YBCO have revealed hysteretic jumps in the 
current- voltage curves at low temperatures, as well as 
non-Gaussian noise fluctuations in the resistance curves 
[16]. As the temperature is increased, the current-voltage 
curves become smooth and non-hysteretic. These results 
have been interpreted as a signature of some form of large 
scale heterogeneities such as disordered fluctuating do- 
mains. Other experiments on cuprate superconductors 
in the non-superconducting region of the phase diagram 
have shown magnetic hysteresis and avalanchelike jumps 
in the magnetization curves [17,18], reminiscent of the 
Barkhausen noise that occurs for domain wall depinning 
in ferromagnets. This magnetic hysteresis vanishes at low 
and high doping and at high temperatures [18]. Trans- 
port experiments in 2DEGs have also produced hysteretic 
current-voltage curves which are believed to indicate the 
formation of pinned clump or bubble phases [19]. 

In this work we study the noise fluctuations and hys- 



teresis in the transport curves for a two-dimensional 
model that exhibits uniform, clump, and disordered 
stripe phases as a function of density. Our model consists 
of interacting particles with a repulsive Coulomb inter- 
action and an additional short range exponential attrac- 
tion. At zero temperature and in the absence of quenched 
disorder, a uniform crystal occurs at low density. As the 
density is increased the system forms a clump phase, then 
a stripe or labyrinth phase followed by an anti-clump 
phase and finally a uniform partially crystalline phase at 
high densities. When quenched disorder is present the 
phases are more disordered and can be pinned in the 
presence of a driving force. 

We find that the heterogeneous clump and stripe 
phases exhibit hysteresis in the transport curves accom- 
panied by a large increase in the noise power with a 1/ f a 
noise spectrum characteristic, where a > 1. In the uni- 
form phase at low and high densities, the hysteresis van- 
ishes and l/f noise appears. As the temperature is in- 
creased in the heterogeneous phases, the hysteresis is lost 
and the transport curves become smooth with an accom- 
panying decrease in the noise fluctuations. We find that 
the hysteresis occurs in a dome-shaped region of the tem- 
perature and density phase diagram. Importantly, these 
properties require both the intrinsic heterogeneity and the 
pinning mechanisms. Our results are consistent with the 
clump and stripe phases forming a glassy or domain glass 
phase. 

Model — We consider a two-dimensional system with 
periodic boundary conditions in the x and y directions. 
We fix the system size to 90 x 90 and vary the particle den- 
sity n by changing the number of particles N from 80 to 
3000. The particle-particle interaction fy = — VU(rij)rij 
for particles separated by ry = rj consists of a long- 
range Coulomb repulsion and a short range exponential 
attraction: 



U(r) — l/r — Bexp(-Kr). 



(1) 



At small and large r the repulsive Coulomb term domi- 
nates. The attractive interaction can be varied using 
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FIG. 1. The particle positions (black dots) after anneal- 
ing to T = for densities (a) n = 0.013, (b) n = 0.06, (c) 
n = 0.15, and (d) n = 0.34. 



the inverse screening length k and the parameter B; in 
this work both B and k are held fixed. In previous 
work we have shown that this model produces crystal, 
clump, stripe, and anti-clump phases as B and the den- 
sity n are varied [20]. Here we also consider the effects 
of quenched disorder, modeled as N p randomly placed 
attractive parabolic pins of strength f p = 2.5 and radius 

r p , giving f p = J2k=i(fp/ r p) r ik®( 1 ~ r tk/r p ), where 6 
is the Heaviside step function, for particle-pin spacing of 
r«fe = r« — ffe- The overdamped equation of motion for 
a single particle i is r\ Vj = f \ = fy + fp + fd + f T - 

Here rj is a phcnomcnological damping term. The driving 
term fd = /^x would come from an applied voltage in the 
case of charged particles. The conduction is proportional 
to the average particle velocity < V > in the direction 
of the applied drive. We increase the drive in small in- 
crements and average over many thousands of time steps 
to avoid transient effects and ensure the velocities have 
reached a steady state. The temperature is modeled as 
Langevin random kicks with the properties < f T (t) >= 
and < f T {t)f T (t') >= 2r)k B T5(t-t'). The system is ini- 
tially prepared in a high temperature molten state and 
annealed down to a lower finite temperature, after which 
the driving force is applied. The parameters we vary in 
this work are the density, driving force, and temperature. 
We have performed a similar series of simulations with 
f p = 1.5 and find the same behavior at shifted values 
of drive. We have also checked the effect of system size 
using a 40 x 40 system, and find the same results shown 




FIG. 2. The velocity V vs driving force f d curves for (a) 
density n = 0.15, T = 0.1; (b) n = 0.06, T = 0.1; (c) n = 0.15, 
T = 1.3; (d) n = 0.34, T = 0.1. 



here but with lower resolution. 

In Fig. 1 we illustrate some of the representative phases 
that occur with increasing density in the presence of the 
quenched disorder. In Fig. 1(a) at density n = 0.013, 
the particles form a uniform phase of single charges. In 
Fig. 1(b) at n = 0.06, the system is comprised of hetero- 
geneous arrangements of disordered clumps. For higher 
density, such as n — 0.15 shown in Fig. 1(c), the system 
forms a disordered labyrinth pattern. At very high den- 
sities the system returns to a uniform phase with consid- 
erable crystalline order as seen in Fig. 1(d) for n = 0.34. 
The heterogeneous phases occur for 0.05 < n < 0.34. 

By examining the velocity-force curves at T = 0, we 
find smooth non-hysteretic curves for densities where uni- 
form phases occur, while in the heterogeneous regions, 
V(fd) shows pronounced hysteresis. In Fig. 2(a) we 
plot the velocity vs applied force curve for a heteroge- 
neous system at n = 0.15 where the equilibrium state 
is a labyrinth phase. Here the curve shows two abrupt 
changes in slope, and V on the decreasing sweep of fd 
is higher than on the increasing sweep. If the drive is 
increased again the same ramp up curve is followed. In 
Fig. 2(b), V(fd) is shown for a system with n = 0.06, 
where the equilibrium system forms a clump phase. Here 
a hysteretic region is also observed. The current-voltage 
curve for n = 0.34, plotted in Fig. 2(d), is smooth and has 
no hysteresis within our resolution. This density n = 0.34 
corresponds to the uniform phase seen in Fig. 1(d). For 
the densities that exhibit hysteresis, as the temperature 
is increased the width of the hysteresis is reduced and it 
vanishes completely at high temperatures. In Fig. 2(c) we 
plot the velocity- force curve for the system in Fig. 2(a) 
for a temperature of T = 1.3 where the hysteresis has 
been lost and the velocity force curve is smooth. 

In order to better characterize where the heterogeneous 
phases occur and where a hysteretic response is observed, 
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FIG. 3. Temperature T vs density n. The shaded region 
indicates where hysteresis (H) is present in the velocity-force 
curve. The nonhysteretic region is indicated by NH. 

we have performed a series of simulations for different 
densities and temperatures. In Fig. 3 we highlight the re- 
gion of temperature and density where the phases exhibit 
a finite hysteresis in the velocity force curves. For high 
and low densities the transport curves are non-hysteretic 
at T = 0. As the density increases, the temperature at 
which the hysteresis disappears increases and reaches a 
maximum at n = 0.15, which corresponds to the density 
at which labyrinth or stripe patterns appear. As the den- 
sity is further increased, the hysteresis width decreases 
and disappears at n « 0.34. There is some asymmetry 
in the hysteretic region, with the hysteresis extending 
further on the higher density side of the peak. This is 
indicative of the importance of collective particle interac- 
tions in producing the hysteresis; there are more particles 
available to interact in the anti-clump state on the high 
density side of the peak than in the clump state on the 
low density side. In recent magnetization experiments 
in doped cuprates, a similar dome-like structure was ob- 
served where hysteresis is present as a function of doping 
and temperature [18] . A direct comparison to the magne- 
tization experiment is difficult since our model does not 
include magnetism; however, hysteresis in magnetic ma- 
terials or superconductors can be modeled in general as 
collections of particles or domains walls interacting with 
quenched disorder and an external driving field, similar 
to the model we are using. It would be interesting to 
measure the hysteresis in the transport characteristics of 
the experimental system as a function of temperature and 
doping to see if a similar dome under which hysteresis is 
present occurs, as found in our simulations. We note that 
hysteresis is not seen in 2D systems with quenched disor- 
der for purely repulsive interactions such as vortices [21]. 
Hysteresis in transport can appear in 3D vortex systems 
when the sample breaks up into two phases and 
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FIG. 4. (a) Time series of the velocity fluctuations for a 
system at T = 0.1 and n = 0.15. (b) The power spectra 
S(v) for the velocity fluctuations. Upper curve: n = 0.15, 
lower curve: n = 0.3. The upper solid line is a power law fit 
with a = 1.6, and the lower solid line is a power law fit with 
a = 0.95. (c) Noise power So for a fixed frequency range vs 
n. (d) Power spectrum exponent a vs n. 

these phases form a heterogeneous labyrinth structure 
[22]. This suggests that hysteresis in transport is a sign 
of large scale heterogeneities. 

We next examine the noise fluctuations at different 
densities and temperatures. We apply a constant drive 
fd and measure the velocity fluctuations SV(t) as a func- 
tion of time, as illustrated in the velocity trace curve of 
Fig. 4(a) for a system with n = 0.15. From the velocity 
fluctuations we can determine the power spectrum using 
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(2) 



The noise power So is defined as the average value of the 
power spectrum over a particular frequency octave. 

In Fig. 4(b) we plot the power spectrum for n = 0.15, 
where the labyrinth phase forms, as well as for n = 0.3 
where the high density uniform phase forms. Both mea- 
surements are performed at T = 0.1. In the labyrinth 
phase a l/f a noise spectrum is observed with a = 1.6, 
while for n = 0.3, a = 1.0. By performing a series of 
simulations, we can plot the noise power and a for varied 
n and T. In Fig. 4(c) we show the noise power So as a 
function of n for fixed T = 0.1. In the hysteretic regions, 
a > 1, while in the uniform regions, a < 1. The noise 
power is maximum near n = 0.15. In Fig. 4(d) we plot 
a, which also shows a maximum in the hysteretic regime. 
The power spectra for the Barkhausen type noise seen in 
the magnetization experiments were not measured. It 
would be interesting to examine the change of the noise 
fluctuations in the hysteretic and non-hysteretic regimes. 
Simulations of systems producing Barkhausen noise give 
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(b) 



experiments in cuprate superconductors [16-18] in which 
charge heterogeneities are believed to occur. 
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FIG. 5. (a) The exponent a vs T obtained from the power 
spectrum for a system with n = 0.15. (b) The corresponding 
noise power 5*0 vs T. 

a = 1.77 [23]. The appearance of a > 1.0 and large 
noise power in experiments and simulations has been in- 
terpreted as evidence of glassy dynamics in 2D charged 
systems [24,25]. In contrast, noninteracting or single 
particle scenarios for noise generation predict a < 1.0 
[26-28]. 

We next examine the noise at a fixed density and in- 
creasing temperature as shown in Fig. 5 for a system with 
n = 0.15. Here the noise power decreases with increasing 
temperature and a also decreases, as seen in Fig. 5(a). 
As the temperature increases, the heterogeneities begin 
to melt and the system becomes more uniform. In gen- 
eral, in the hysteric regions of the T-n phase diagram, 
a > 1. In Fig. 5(a) the hysteresis disappears in the 
velocity force curves for T > 1.0. This corresponds to 
the temperature regime where the value of a drops to 
a < 1. The noise power So decreases exponentially with 
temperature, as shown in Fig. 5(b). Single particle hop- 
ping models often predict an increase in the noise power 
with temperature rather than the decrease observed here 
[28] . This supports the importance of collective particle 
interactions in determining the transport properties and 
hysteresis. 

In summary, we have examined hysteresis and noise 
in a model system of heterogeneous charge ordering in 
the presence of quenched disorder. We find that at low 
and high densities, the system forms non- heterogeneous 
phases that do not exhibit hysteresis in the transport 
curves. At intermediate densities, however, the system 
forms heterogeneous clump and stripe phases which ex- 
hibit hysteresis in the transport curves. As a function 
of temperature and doping, we find a dome region where 
the hysteresis occurs. This prediction can be tested in, 
e.g., transport measurements in high temperature su- 
perconductors for varied dopings. We also examine the 
noise properties and find that in the heterogeneous and 
hysteretic regions, the noise power is substantially en- 
hanced and shows a l/f a noise spectrum characteristic 
with a > 1 . We predict a peak in the noise power and a 
as a function of particle density. Additionally both the 
noise power and the exponent a in the heterogeneous re- 
gions decrease as a function of temperature. Our results 
are consistent with recent transport and magnetization 
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